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Abstract 



I summarize our self-contained determinations of the lowest order hadronic contributions (l], |^] to the 
anomalous magnetic moments a MjT of the muon and tau leptons, the running QED coupling a(Mz) and 
the muonium hyperfine splitting v. Using an average estimate of the light-by light scattering contribu- 
tion: a^LL) = 85(18) x lO" 11 , we deduce: a^ M = 116 591 861(78) x HT 11 , af M = 117 759(7) x lO" 8 , 
giving: a^ M -a^ p = 162(170) x 10 -11 . We also obtain: a~ 1 {M z ) = 128.926(25) and the Fermi energy 
splitting: Vp M = 4 459 031 783(229) Hz. Lower bounds on some new physics are given, while vf. M 
leads e.g. to m Pl jm e = 206.768 276(11) in remarkable agreement with the data. 

o \ 

1 Introduction 

Oh! 

QED is at present the gauge theory where perturbative calculations are the most precise known 
today. Therefore, accurate measurements of QED processes are expected to give strong constraints 
on different electroweak models and to reveal some eventual deviations from the standard model (SM) 
predictions ^. In the following, I will discuss the effects of the hadronic and QCD contributions to three 
classical QED processes which are: the anomalous magnetic moment of the muon and tau leptons, the 
running QED coupling a(Mz) and the muonium hyperfine splitting v. These hadronic contributions 
are one of the main sources of uncertainties into these processes. Using a dispersion relation, it is 
remarkable to notice that the different lowest order hadronic contributions for these three processes 
can be expressed in a closed form as a convolution of the e + e~ — > hadrons cross-section an{i) with a 
QED kernel function K{t) which depends on each observable: 

Chad = A / dt R o(t) <rnQt) , (1) 



47T 3 

where: 



4/jiz 



Chad = a; ih ad , Aa ha d x 10 5 or Av had . (2) 



x Talk presented at the 1st High-Energy Physics Madagascar International Conference Series (HEP-MAD'01), 27th 
sept-5th Oct. 2001, Antananarivo (to be published by World Scientific, Singapore). 

2 For general discussions on astroparticle physics and supersymmetric models, see e.g. ji| ^. 
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For the anomalous magnetic moment a/ h a( j, K ai (t > 0) is the well-known kernel function H: 

Ka,(t) = I dX x 2 + (t/mf)(l-x) ' (3) 
where mi is the lepton mass. It behaves for large t as: 

2 

771 

# 0! (*> m 2)~_i- . (4) 
For the QED running coupling Aa had x 10 5 , the kernel is (see e.g. M): 

n 

where a _1 (0) = 137.036 and Mz = 91.3 GeV. It behaves for large t like a constant. 
• For the muonium hyperfine splitting Ai/had) the kernel function is (see e.g 

K v = -p, 

where: 



(x^ + 2) In j—f- - ^ + | ) ln% 



(0) 




m e £ /I 

and we take (for the moment) for a closed comparison with R| |^], the value of the Fermi energy 
splitting: 

u F = 445 903 192 0.(511)(34) Hz . (8) 

It behaves for large t as: 

/ Tn 2 \ /n + it;\ 

(9) 

The different asymptotic behaviours of these kernel functions will influence on the relative 
weights of different regions contributions in the evaluation of the above integrals. 

2 Input and Numerical Strategy 

The different data input and QCD parametrizations of the cross-section an(t) have been discussed 
in details in Q] (herereferred as SN1) and corresponding discussions will not be repeated here. The 
sources of these data are quoted in the last column of Table 1 from SN1 and Table 2 from [| (her- 
ereferred as SN2) are classified according to the estimate in different regions. We shall only sketched 
briefly the numerical strategy here: 

• Our result from the 1 = 1 isovector channel below 3 GeV 2 is the mean value of the one using 
r-decay and e + e~ data. In both cases, we have used standard trapezoidal rules and/or least 
square fits of the data in order to avoid theoretical model dependence parametrization of the 
pion form factor. Correlations among different data have been taken in the compilations of Q 
used in this paper. In the region (0.6 — 0.8) GeV 2 around the u-p mixing, we use in both cases 
e + e~ data in order to take properly the SU(2)p mixing. The SU(2) breaking in the remaining 
regions are taken into account by making the average of the two results from r-decay and e + e~ 
and by adding into the errors the distance between this mean central value with the one from 
each data. 



3 In the next section, we shall extract this value from the analysis. 
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• For the 1 = isoscalar channel below 3 GeV 2 , we use the contributions of the resonances 
w(782) and </>(1020) using narrow width approximation (NWA). We add to these contributions, 
the sum of the exclusive channels from 0.66 to 1.93 GeV 2 . Above 1.93 GeV 2 , we include the 
contributions of the w(1.42), cj(1.65) and </>(1.68) using a Breit-Wigner form of the resonances. 

• For the heavy quarkonia, we include the contributions of known J/tp (IS to 4.415) and T (IS 
to 11.02) families and use a NWA. We have added the effect of the it bound state using the 
leptonic width of (12.5 ± 1.5) keV given in ||. 

• Away from thresholds, we use perturbative QCD plus negligible quark and gluon condensate 
contributions, which is expected to give a good parametrization of the cross-section. These 
different expressions are given in SN1. However, as the relative role of the QCD continuum is 
important in the estimate of Aahad, we have added, to the usual Schwinger interpolating factor 
at order a s for describing the heavy quark spectral function, the known a^m,Q/t corrections 
given in SN1. However, in the region we are working, these corrections are tiny. 

• On the Z-mass, the integral for Aahad has a pole, such that this contribution has been separated 
in this case from the QCD continuum. Its value comes from the Cauchy principal value of the 
integral. 

3 Lowest order hadronic contributions 

3.1 Muon and tau anomalies 

We show in Table 1 the details of the different hadronic contributions from each channels and from 
different energy regions for the muon and tau anomalies. Taking the average of the results in Table 1 
and adding further systematics due to an eventual deviation from the CVC assumption and from the 
choice of the QCD continuum threshold for the light flavours, one deduce the final estimate from an 
average of the e + e~ and r-decay data 

a* ad (l.o) = 7020.6(75.6) x 10~ n , a h T ad {l.o) = 353.6(4.0) x 10~ 8 , (10) 

• The main error (80% when added quadratically) in our previous determinations comes from the 
p-meson region below 0.8 GeV 2 . Hopefully, improved measurements of this region are feasible 
in the near future. 

• The second source of errors comes from the region around M T for the inclusive r-decay and 
between 1 GeV to M T for the e + e~ data. These errors are about half of the one from the 
region below 0.8 GeV 2 in most different determinations. They can be reduced by improved 
measurements of inclusive r-decay near M T (1 = 1) and by improving the measurements of the 
odd multi-pions and KK, KKn, ... channels in the 1 = channels from e + e~ data. 

• The contributions of the whole region above 3 GeV 2 induce much smaller errors (7% of the 
total). There is a quite good consensus between different determinations in this energy region. 

• These predictions agree within the errors with previous predictions quoted in SN1 Q and recent 
estimates given in Q-Q. 

3.2 Running QED coupling 

Using the same data as for the anomalous magnetic moment, one can deduce from Table 2: 

Aa had = 2763.4(16.5) x 10~ 5 . (11) 

Also a detailed comparison of each region of energy with the most recent work of (9| shows the same 
features (agreement and slight difference) like in the case of in SN1 due to the slight difference 
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in the parametrization of the data and spectral function. However, the final results are comparable. 
Finally, one can remark that due to the high-energy constant behaviour of the QED kernel function 
in this case, the low-energy region is no longer dominating. For a^, the contribution of the p-meson 
below 1 GeV is 68% of the total contribution, while the sum of the QCD continuum is only 7.4% (see 
e.g. SN1). Here the situation is almost reversed: the contribution of the p-meson below 1 GeV is only 
2%, while the sum of the QCD-continuum is 73.6%. 



Table 1: 



Determinations of a^ ad 



columns) and averaged e + e~ 



(l.o) using combined 
data (3rd column). 



e and inclusive r decay data (2nd and 4th 



Region in GeV 2 



afT d (l.o) x 10 1 



„had 



l.o) x 10 s 



Data input 



Light Isovector 

4m 2 -> 0.8 
0.8 -> 2.1 
2.1 -> 3. 
Total Light 1=1 
Light Isoscalar 
Below 1.93 



r+e + e 

4794.6 ± 60.7 
494.9 ± 15.8 
202.0 ±29.7 
5491.5 ± 69.4 



4730.2 ± 99.9 
565.0 ± 54.0 
175.9 ± 16.0 
5471.1 ± 114.7 



r+e + e 

165.8 ± 1.5 
28.7 ± 1.1 
17.0 ±2.6 
211.5 ± 3.2 





387.5 ± 13 


387.5 ± 13 


15.3 ±0.5 


4> 


393.3 ± 9.9 


393.3 ±9.9 


21.0 ±0.5 


0.66 -» 1.93 


79.3 ± 18.7 


79.3 ± 18.7 


4.3 ± 1.1 


From 1.93 to 3 








w(1.42), w(1.65) 


31.3 ±6.8 


31.3 ±6.8 


2.6 ±0.7 


0(1.68) 


42.4 ± 18.2 


42.4 ± 18.2 


3.8 ± 1.3 


Total Light 1=0 


933.8 ± 31.5 


933.8 ± 31.5 


47.0 ± 2.0 


Heavy Isoscalar 








J/ip(lS -» 4.415) 


87.0 ± 4.7 


87.0 ±4.7 


13.08 ±0.69 


T(1S -> 11.020) 


0.95 ± 0.04 


0.95 ± 0.04 


0.23 ±0.01 


Total Heavy 1=0 


88.0±4.7 


88.0±4.7 


13.3 ± 0.7 


QCD continuum 








3. -> (4.57) 2 
(4.57) 2 -> (11.27) 2 
(11. 27) 2 -» 4M t 2 
4M t 2 -> 00 


407.0 ± 2.3 


407.0 ± 2.3 


49.4 ±0.3 


95.3 ±0.5 


95.3 ±0.5 


27.3 ±0.1 


20.5 ±0.1 


20.5 ±0.1 


5.87 ±0.01 


w 0. 


w 0. 


« 0. 


Total QCD Cont. 


522.8 ±2.4 


522.8 ± 2.4 


82.6 ± 0.3 




NWA [14 1 
NWA jjlj 

exclusive |l£ 



BW 
BW 



NWA 
NWA 



[14| 
[14| 



(u,d,s) 
(u, d, s, c) 
(u, d, s, c, b) 
(u, d, s, c, b, t) 



7036.1(76.4) 



7015.7(119.1) 



354.4(3.; 



For this reason, improvement due to the new Novosibirsk e + e~ data in the low-energy region is not 
significant, as we have explicitly checked. At present, new BES data |^| in the J/tp region are also available, 
which can be alternatively used. Below the J/tp resonances, the BES data are in excellent agreement with 
the QCD parametrization to order used here for 3 flavours, justifying the accuracy of your input. Above 
the J/ip resonances, the parametrization used here (sum of narrow resonances ±QCD continuum away from 
thresholds) can also be compared with these data. On can notice that, in the resonance regions, the BES 
data are more accurate than previous ones, which may indicate that our quoted errors in Table 1 for the 
J jip family contributions are overestimated. In addition, the threshold of the QCD continuum which we have 
taken above the J/ip resonances, matches quite well with the one indicated by the BES data. Our estimate of 
Aahad is compared with recent predictions respectively in Fig. 1 from SN2 0] (references to the authors are 
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20|) where one can notice a very good agreement. 



Table 2: Lowest order determinations of Aahad x 10 5 and Ai/h a d [Hz] using combined e + e and 
inclusive r decay data (2nd and 4th columns) and averaged e + e~ data (3rd and 5th columns). 



Region in GeV 2 



Aa had x 10 5 



A^ had [Hz] 



Data 



Light Isovector 

4m 2 -> 0.8 
0.8 -> 2.1 
2.1 3. 

Total Light 1=1 
Average 

Light Isoscalar 

Below 1.93 

<t> 

0.66 -» 1.93 
From 1.93 to 3 
w(1.42), w(1.65) 
0(1.68) 

Heavy Isoscalar 

J/V>(15 -> 4.415) 
T(LS-> 11.020) 
T(349) 

Total Heavy 1=0 
QCD continuum 

3. -> 4.57 2 

4.57 2 -> 11. 27 2 

11.27 2 -> (M z -3 GeV) 2 

(M z + 3 GeV) 2 -> 4M t 2 

Z-pole 

4M t 2 -> oo 

Toia/ QCD Coni. 



r+e e 

314.5 ± 2.3 

77.2 ±3.4 

62.3 ±9.2 
454.0 ± 70.6 



e e 

302.7 ±7.1 
82.0 ±5.4 
53.6 ± 4.9 



r+e e 

152.9 ± 1.8 
12.1 ± 0.5 
7.8 ± 1.2 



e e 

148.4 ±3.1 
16.9 ± 1.9 
6.7 ±0.6 



.2 ±10.2 172.8 ±2.2 112.1 ±3.1 



446.1 ± 10.4 ± 7.9 



31.5 ± 1.1 
52.3 ± 1.2 

11.6 ±3.0 

9.4 ± 1.4 

14.6 ±4.6 

119.0 ±5.9 

116.3 ±6.2 

12.7 ±0.5 
-(0.1 ±0.0) 
128.9 ±6.2 

330.1 ± 1.0 
503.0 ± 1.0 
2025.7 ± 2.0 

-(794.6 ±0.6) 

29.2 ± .5 
-(24.0 ±0.1) 
2069.4 ± 5.2 



112.5 ±3.0 ±0.3 



12.7 ±0.4 

13.7 ±0.3 
2.7 ±0.7 

1.2 ±0.2 
1.7 ±.5 

32.1 ± 1.0 

4.0 ±0.2 
0.1 ±0.0 

w 

4.1 ± .2 

17.5 ±.1 
5.0 ±.1 

1.3 ±0.0 

w 
w 
« 

23.8 ± 1.4 




NWA [14 1 
NWA jw| 
excl 



BW 
BW 



[16, 14 
[16, 17, 



NWA 
NWA 
NWA 



1141 
1141 
[14, 



(u,d,s) 
(u, d, s, c) 
(u, d, s, c, b) 

ppal value || 
(u, d, s, c, 6, i) 



Final value 



2763.4 ± 16.5 



232.5 ±3.2 



3.3 The muonium hyperfine splitting 

Our final result for A^h a d comes from Table 2 of SN2 

A^ ha d = (232.5 ± 3.2) Hz 



(12) 



and is shown in Table 3 in comparison with other determinations, where there is an excellent agreement with 
the most recent determination M . Here, due to the (In t) ft behaviour of the kernel function, the contribution 
of the low-energy region is dominant. However, the p-meson region contribution below 1 GeV is 47% compared 
with 68% in the case of a^, while the QCD continuum is about 10% compared to 7.4% for a M . The accuracy 
of our result is mainly due to the use of the T-decay data, explaining the similar accuracy of our final result 
with the one in [Q| using new Novosibirsk data. The agreement with Q can be understood from the agreement 
of the averaged correlated e + e" and r-decay data compiled in J§| with the new Novosibirsk data used in Q. 
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10 5 Aa h 



ADH 98 
DH 98 
BP 01 

MOR 01 

J 01 

YT01 

SN01 
(this work) 



a -1 (M z ) 



2700 2800 



128.85 128.90 128.95 129.00 



Figure 1: Recent determinations of A«had and a 1 (Mz). The dashed vertical line is the mean central 
value. 



However, we differ with DH98 @ in the treatment of the QCD contribution ^ For light quarks, QCD is 
applied in the region where non-perturbative contributions are inessential. For heavy quarks, perturbative 
QCD is applied far from heavy quark thresholds, where it can be unambiguously used. 



Table 3: Recent determinations of Ai^ 



.. i 



Authors 



Afhad [Hz] 



FKM 99 Jp 
CEK 01 ff 

sn oi fli 



240 ±7 
233 ±3 
232.5 ±3.2 



4 Theory of the muon and tau anomalies and new physics 
4.1 QED and EW contributions to a M 



The QED up to 8th order and EW including two-loop corrections n are: 



QED 



116 584 705.7(2.9) x 10" 



151(4) x 10" 



(13) 



4.2 Higher order hadronic contributions to a. 



Higher order hadronic contributions have been discussed first in j22j. They can be divided into two classes 
The first one involves the vacuum polarization and can be related to the measured e + e~ — ► hadrons total 
cross-section, similar to the lowest order contribution. After rescaling the result in bA M, one obtains pi]: 



had / 7 \ 

a„ {n.o.jv.p 



-101.2(6.1) x 10" 



(14) 



4 For more details, sec . 

5 References to original works can be found in [M pi. 
6 For more details of the following discussions, see [b 
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The second class is the light-by-light scattering diagram. Contrary to the case of vacuum polarizations, this 
contribution is not yet fully related to a direct measurable quantity. In order to estimate this contribution, 
one has to introduce some theoretical models. The ones used at present are based on chiral perturbation J2f| 
and ENJL model To both are added vector meson dominance and phenomenological parametrization 

of the pion form factors. The different contributions are summarized in Table ^j, where the first two come 
from the boson and quark (constituent) loops, while the last one is due to meson pole exchanges. The first 
two contributions are quite sensitive to the effects of rho-meson attached at the three off-shell photon legs 
which reduce the contributions by about one order of magnitude (!). The third one with pseudoscalar meson 
exhanges (anomaly) gives so far the most important contribution. There is a complete agreement between 
the two model estimates (after correcting the sign of the pseudoscalar and axial- vector contributions J27j), 
which may indirectly indicate that the results obtained are model-independent]^. Neverthless, there are still 
some reamining subtle issues to be understood (is the inclusion of a quark loop a double counting ?, why the 
inclusion of the rho-meson decreases drastically the quark and pion loop contributions ? is a single meson 
dominance justified?..). The results in Refs. and after correcting the sign of the pseudoscalar and 
axial- vector contributions as suggested in p7[|, are given in the table: 



Table 4: a£ ad (h.o) LL x 10 11 



Type of diagrams Ref. [25| Ref. 



tt~ loop -4.5(8.1) -19(13) 

quark loop 9.7(11) 21(3) 

7T°, Vl rj' poles 82.7(6.4) 85(13) 

axial-vector pole 1.74 2.5(1.0) 

scalar pole *> -6.8(2.0) 

Total 82.8(15.2) 82.7(18.8) 

*^ Wc have added here the result from . 

while a naive constituent quark model gives || using the result of [^] : 

a^ ad (h.o) LL = 92(20) x 10" 11 . (15) 

Due to the unknown real value of the virtual photon momenta entering into the calculation, it can happen that 
none of the previous approaches describes accurately the LL effect |^. Therefore, for a conservative estimate, 
we take an arithmetical average of the three determinations: 

a% ad (h.o) L L = 84.7(18.0) x 10" 11 . (16) 

One can notice, for the muon, a strong cancellation between the higher order vacuum polarisation and the 
light by light scattering contributions. 

4.3 The total theoretical contributions 

Summing up different contributions, the present theoretical status in the standard model is Jl|: 

al M = 116 584 840.2(19.6) x lO^ 11 + a% ad (l.o) 

= 116 591 860.8(78.1) x 10" 11 , (17) 

7 See however 

8 I thank Eduardo de Rafael and Paco Yndurain for some clarifying communications on this point. 
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where a^ ad (Lo) is the lowest order hadronic contributions evaluated in SN1 § (see Eq. (0)). Comparing this 
SM prediction with the measured value Jlo| ]: 

= 116 592 023(151) x 10" 11 , (18) 

we deduce: 

a" em = al xp - al M = 162(170) x 10~ n . (19) 

If the future data will be accurate by ±40 x 10 -11 , while the theoretical errors are almost unchanged due 
to different limitations discussed previously, then the error in a^ ew will be reduced by a factor 2, giving a 
chance to detect a 2a deviation from SM if any. Combined with the mean value of existing determinations of 
a ti (l.o) had given in SN1 jl[, which gives: a™"" = 175(170) x 10" 11 , one can deduce the conservative range: 

- 56 < a n ™ x 10 11 < 393 (90% CL) . (20) 

4.4 Bounds on some new physics from a M 

This result gives, for a supersymmetric model with degenerate sparticle mass, a lower bound of about 113 
GeV |j, while the compositeness and the leptoquark scales are constrained to be above 1 TeV. Bound on the 
sparticle mass is comparable with present experimental bound from LEP data. The one of the leptoquarks is 
much larger than the present lower bounds of about (200 ~ 300) GeV from direct search experiments at HERA 
and Tevatron. We expect that these different bounds will be improved in the near future both from accurate 
measurements of a M and of e + e~ data necessary for reducing the theoretical errors in the determinations of 
the hadronic contributions, being the major source of the theoretical uncertainties. 

4.5 Tau anomaly 

In the same way, the higher order hadronic contributions read [jl], |3l] | 

a h ^ ad {h.o)vp = 7.6(0.2) x 10~ 8 , a^ d (h.o) LL = 23.9(5.1) x 10" 8 , (21) 
which, in the r case, the two effect add each others. Adding the other contributions, one obtains M; 

af M = 117 759.1(6.7) x 10 -8 . (22) 
This value can be compared with the present (inaccurate) experimental one |B2|: 



exp 



0.004 ± 0.027 ± 0.023 , (23) 



which, we wish, will be improved in the near future. 



5 The QED running coupling a(Mz) 



To the lowest order hadronic contribution in Eq. (|ll[), we add the radiative corrections taken by adding the 
effects of the radiative modes 7r°7, 777, ■K + -K~ r y, .... We estimate such effects to be: 

Aa had = (6.4 ± 2.7) x KT 5 (24) 

by taking the largest range spanned by the two estimates in Q and [g]. Using the QED contribution to 
three-loops |J: 

Aqqed = 3149.7687 x 10~ 5 , (25) 
and the Renormalization Group Evolution of the QED coupling: 

a -1 (M z ) =a _1 (0) 1- Aqqed - Aa had , (26) 

one obtains the final estimate: 

oT 



_1 (M Z ) = 128.926(25) , (27) 

r ith recent existing determinations. One can notice 
>ns [|, [], po|| . We expect that with this new impro 
of a(Mz), present lower bound of 114 GeV from LEP data on the Higgs mass can be improved. 



which we show in Fig 1 for a comparison with recent existing determinations. One can notice an improved 
accuracy of the different recent determinations j^, [|, [], po|| . We expect that with this new improved estimate 



9 For more recent detailed discussions, see e.g. |i[ Q and references quoted there. 
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6 Muonium hyperfine splitting 

Adding to this result in Eg. ( |l^ ) from SN2 || , the QED contribution up to fourth order, the lowest order 
electroweak contribution Q, and an estimate of the higher order weak and hadronic contributions [Q: 

A^qed = 4 270 819(220) Hz ,Ai/ weak (Lo) = -^|m e m M v F ~ -65 Hz , 

|Ai/ wea k(MI 0.7 Hz , Ai^ad(M ~ 7(2) Hz , (28) 

one obtains the Standard Model (SM) prediction: 

V S M = Vp + A^QED + Ai/ wcak + A^ had + Av had (h.o) . (29) 

If one uses the relation: 

= " (t) dW^) 3 : " = ' (30) 

and Z = 1 for muonium, a _1 (0) = 137.035 999 58(52) @, cR^ =3 289 841 960 368(25) kHz @, one would 
obtain: 

u SM = 4 463 302 913(511) (34)220) Hz , (31) 

where the two first errors are due to the one of the Fermi splitting energy. The first largest one being induced 
by the one of the ratio of the magnetic moments. The third error is due to the 4th order QED contribution 
where, one should notice that, unlike the case of a M , the dominant errors come from the QED calculation 
which should then be improved. This prediction can be compared with the precise data pfj]: 

z, oxp = 4 463 302 776(51) Hz . (32) 

Therefore, at present, we find, it is more informative to extract the Fermi splitting energy vf from a comparison 
of the Standard Model (SM) prediction with the experimental value of v. Noting that vf enters as an overall 
factor in the theoretical contributions, one can rescale the previous values and predict the ratio: 

— = 1.000 957 83(5) . (33) 

VF 

Combining this result with the previous experimental value of v, one can deduce the SM prediction: 

uf M = 4 459 031 783(226) Hz , (34) 

where the error is dominated here by the QED contribution at fourth order. However, this result is a factor 
two more precise than the determination in jjj given in Eq. (|^), where the main error in Eq. ([]) comes from 
the input values of the magnetic moment ratios. Using this result in Eq. (^) into the expression: 

vf=Pf(— ) t— —, rg(l + ev) , (35) 

where pF is defined in Eq. (|3c|), and a M = 1.165 920 3(15) x 10 -3 (n]], one can extract a value of the ratio of 
the muon over the electron mass: 

^ = 206.768 276(11) , (36) 

m e 

to be compared with the PDG value 206.768 266(13) using the masses in MeV units, and with the one from 
0: 206.768 276(24). In Ref. an accuracy two times better than the present result has been also obtained. 
However, in that case, the errors in the QED contributions may have been underestimated. After inserting 
the previous value of m e /m M into the alternative (equivalent) relation given in Eq. (po[), one can deduce the 
ratio of magnetic moments: 

^- = 4.841 970 47(25) x 10" 3 , (37) 
Ms 

compared to the one obtained from the PDG values of p u l Pv and ^ p //i fl [Q: fi^/fXg = 4.841 970 87(14) x 
10 -3 . In both applications, the results in Eqs. (^) and ( p7[ ) are in excellent agreement with the PDG values. 
These remarkable agreements can give strong constraints to some contributions beyond the Standard Model 
and are interesting to be explored. 
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7 Conclusions 



We have evaluated the lowest order hadronic and QCD contributions a^ ad (l.o), A«h a d an d A^had respectively 
to the anomalous magnetic moment, QED running coupling and to the muonium hyperfine splitting. Our 
self-contained results given in Eqs. @, (0) and @, obtained within the same strategy and data input, are 
in excellent agreement with existing determinations and are quite accurate. We have also revised the estimate 
of the light by light scattering contributions to a^ tT , and have explored some phenomenological consequences 
of these results. 
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